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Abstract 

We present a general derivation the three-body spectator (Gross) equations and the correspond- 
ing electromagnetic currents. As in previous paper on two-body systems, the wave equations and 
currents are derived from those for Bethe-Salpeter equation with the help of algebraic method using 
a concise matrix notation. The three-body interactions and currents introduced by the transition 
to the spectator approach are isolated and the matrix elements of the e.m. current are presented 
in detail for system of three indistinguishable particles, namely for elastic scattering and for two 
and three body break-up. The general expressions are reduced to the one-boson-exchange approx- 
imation to make contact with previous work. The method is general in that it does not rely on 
introduction of the electromagnetic interaction with the help of the minimal replacement. It would 
therefore work also for other external fields. 

PACS numbers: 11. lOSt, 13.40.-1,21.45. +v,24.10.Jv 
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I. INTRODUCTION 



Exposing few nucleoli systems to medium energy (~ lGeV) electromagnetic probes re- 
quires a relativistic description of initial and final states as well as of the reaction mechanism. 

fl 

The relativistic effects are usually included in perturbative way but even at leading order 
one has to include large number of terms into NN potential and electromagnetic current, 
while the convergence of the series in v/c ~ p/m becomes questionable. The conceptually 
cleaner alternative is to start from an explicitly covariant formulation. A number of covari- 
ant techniques has been introduced and studied extensively for the case of the simplest two 
nucleon system, in particular in calculations of the elastic e.m. deuteron form factors |2| and 
in deuteron electro- and photodisintegration 



a a. 



In the present paper we adopt the spectator (or Gross) quasipotential approach 
It reduces four dimensional integrals over loop momenta to their more treatable three di- 
mensional form in explicitly covariant way by picking up the poles of nucleon propagators. 
This simplifies an analytic structure of dynamic equations and leads to closer approximation 
to full result compared to the Bethe-Salpeter ladder approximation. A set of realistic one- 
boson-exchange potentials were constructed within the formalism, giving a fair fit to NN 
scattering data and deuteron properties f|. It was found that the model gives very good 
description of the elastic deuteron form factors 

aa 

. Three nucleon spectator equation with 
those realistic potentials was solved by Stadler and Gross [8] and the normalization condition 
for the 3N bound state amplitudes was derived 0. The first calculation of electromagnetic 
break-up of 3N bound state are in progress. 

In order to calculate this process it is necessary to determine the correct effective elec- 
tromagnetic current operator to be used with the three-body spectator equation. So far 
this problem has been approached in two ways. The first of these is the elegant method of 



Kvinikhidze and Blankleider 







11 



121 ] that gauges the integral equations which describe 



the strongly interacting system in order to obtain the electromagnetic current matrix ele- 
ments for various processes. The second approach is that of Gross, Stadler and Pena [l^ 
which uses an approach based on a diagrammatic description of the three-body system and 
relies heavily on the gauging approach. These two approaches lead to identical results once 
differences in organization are accounted for. In both of these papers certain assumptions 
are inherent in the results that are presented Q|. In both cases, no explicit three-body 
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forces are included and the results implicitly assume that the interaction kernels are given 
by the one-boson-exchange approximation. 

In a previous paper 14] we pursued a different strategy to obtain a general derivation of 
the e.m. currents for the two-body Gross equation. The starting point was the organization 
of the corresponding Feynman amplitudes into 2N Bethe-Salpeter equation and factorization 
of five point Green function (with the external e.m. field). The spectator equations were then 
obtained by separating the propagator of one of the two particles into a piece containing only 
the on-shell positive-energy contribution and a remainder. The Bethe-Salpeter amplitudes 
are then rearranged into integral equations which contain only the on-shell propagator and 
a quasipotential equation which sums all of the off-shell contributions. The electromagnetic 
current consistent with this rearrangement was derived with the help of operator algebra, 
which defines the products of singular operators and represents an useful shortcut alter- 
native to pole analysis of the Feynman diagrams. We have shown that this definition is 
consistent with gauge invariance when a necessary truncation of the strong kernel is applied 
in any finite order of quasipotential decomposition and in a number of meson exchanges 
involved. The main aim of the present paper is to extend this treatment to 3N spectator 
equation and to derive general expressions for electromagnetic matrix elements needed for 
the future application to actual calculation of the electromagnetic properties of relativistic 
three nucleon systems. 

For convenience, Sec. II contains a brief review of the application of this approach to the 
two-body system. 

The starting point of the analogous development for three-body systems is a description 
of the three-body Bethe-Salpeter amplitudes given by four-dimensional integral equations 
that iterate kernels consisting of all three-body irreducible contributions to the six-point 
function. The currents consistent with these kernels can be obtained either by attaching a 
photon to all internal particle propagators in a diagrammatic representation the kernels, or 
by use of the gauging method [12]. Section III contains a review of the three-body Bethe- 
Salpeter equations and current including contributions from two- and three-body kernels 
and currents. Since the three-body spectator equations are defined only for three identical 
particles a discussion of the three-body equations for identical particles is given in Sec. IIIC. 
Before using our approach to obtain the spectator equations, it is convenient to use the 
symmetry of the various Faddeev amplitudes to limit the number of channels that explicitly 
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appear in the equations. We found that it is most convenient to do this in the context of 
a matrix representation of the channel space. Appendix A describes the matrix method for 
three distinguishable particles and Appendix B describes the case of identical particles and 
show how to reduce the problem to two channels. 

Section IV applies the operator technique to the two-dimensional matrix representation 
of the Bethe-Salpeter equation to obtain the corresponding spectator equations for the t- 
matrix and wave functions, and obtains an expression for the corresponding effective current 
operator. Some care is taken to separate two-body quasipotential and scattering matrix 
contributions from three-body contributions that are generated by the transformation from 
the Bethe-Salpeter to the spectator equations. These contributions can occur even if there 
are no explicit three-body Bethe-Salpeter interactions. 

In Sec. V we show how the general expressions obtained in Sec. IV can be reduced to the 
case of the one-boson-exchange approximation for the Bethe-Salpeter kernels and currents. 
This gives results identical to those of Gross, Stadler and Pena; and, by extension, to those 
of Kvinikhidze and Blankleider. 



II. REVIEW OF THE TWO-BODY EQUATIONS 

The motivation for the arguments that we present below for constructing the three- 
body electromagnetic current matrix elements for the three-body spectator equation is best 
illustrated by providing an outline of the arguments that appeared in our previous paper 
|t3 | for obtaining the equivalent expressions for the two-body spectator equation. 

First consider the two-body Bethe-Salpeter equation for two particles interacting by the 
exchange of a boson. Examination of all of the diagrams contributing to the two-body t- 
matrix (the four-point vertex function) allows to show that the sum of all contributions to 
the t-matrix can be represented by an integral equation of the form 

M = V -VG BS M , (1) 

where the kernel of the equation V is the sum of all two-body irreducible diagrams, and where 
Gbs — —iG\G<i is the free two-body Bethe-Salpeter propagator. The formalism reviewed in 
this section and further developed in this paper for three particle systems is restricted to the 
energies below particle production threshold. That is, the propagation of the constituents is 
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approximated by free Green's functions with the physical mass and the phenomenologically 
dressed interaction vertices are assumed, i.e., explicit loops dressing two- and three-point 
Green's functions and their renormalization are not considered. The two-body interacting 
propagator (four-point propagator) can then be written as 

Q = G B s - GbsMGbs = G B s - GbsVQ = G B s - QVG B s ■ (2) 

This propagator contains all of the information about the two-body system and the bound 
and scattering states wave functions can be obtained from the singularities of this propagator. 
By looking for sub-threshold poles in the propagator the bound state Bethe-Salpeter wave 
function can be determined to be 

W) = g bs |r> , (3) 

where |T) is the bound-state Bethe-Salpeter vertex function which satisfies the equation 

ID = -vGbs |r> . (4) 

The pole in the propagator associated with the two incoming particles being on mass shell 
gives the scattering wave function with outgoing wave boundary conditions 

= (1 - G B sM)\p 1 ,s 1 ;p 2 ,s 2 ) (5) 

and the corresponding pole for the outgoing particles gives the conjugate scattering wave 
function with incoming wave boundary conditions 

(^\ = (p 1 ,s 1 ;p 2 ,s 2 \(l-MG B s) ■ (6) 
All of the wave functions satisfy the wave equation 

g- 1 \^) = o. (7) 

The relationship of the two-body interacting propagator to the wave functions implies 
that we can obtain all of the necessary information concerning the electromagnetic currents 
for this system in one-photon-exchange approximation by examining all Feynman diagrams 
contributing to a five-point propagator with one photon leg and four interacting particle 
legs. This five-point propagator can be written as 

gi* = -gj»g (8) 
5 



by defining the current operator 

JM = ii7 M G -l + ij£ Gl l + jW* f (9) 

where Jf is the one-body current operator for particle % and ^ s the two-body current 
composed of all two-body irreducible diagrams involving the absorption of a photon. Any 
current matrix element can then be obtained by picking out the poles in the initial and final 
interacting propagators corresponding to either bound or scattering states. This current 
operator satisfies the Ward-Takahashi identity 

q,J t "=[e 1 (q) + e 2 (q),g- 1 ] , (10) 

where is the photon momentum, e t (q) is the product of the charge e t for particle t (which 
might be an operator in isospin space) and a four-momentum shift operator defined such 
that 

(p't\e t (q)\Pt}= e t (27r)W (p' t -p t -q). (11) 

This along with the wave equation guarantees that the current will be conserved for all 
possible two-body matrix elements. 

The spectator or Gross equation (for distinguishable particles) can be obtained from 
the Bethe-Salpeter equation by substituting a two-body noninter acting propagator gi with 
particle 1 on mass shell for the Bethe-Salpeter propagator Gbs- The t-matrix will remain 
unchanged if it is written as a pair of coupled equations 

M — U — Ug\M (12) 

and 

U = V - VA gi U , (13) 

where the new kernel U is called the quasipotential and Agi = Gbs ~ 9i- 

In ^| it was shown that there exist a set of relationships for the spectator scattering 
matrix, propagators, wave functions and currents which are of the same form as those given 
above for the Bethe-Salpeter equation and that an effective current operator can be derived, 
satisfying a Ward-Takahashi identity which leads to current conservation when used with 
the spectator wave functions. For identical particles the formalism can be represented in a 
compact and transparent matrix notation. 
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We will now show that a similar organization can be obtained for the three-body Bethe- 
Salpeter and Gross equations. We will do this assuming that there may be explicit three- 
body interactions in addition to two-body interactions and will obtain expressions for the 
three-body Gross equations and effective currents which can be expanded to arbitrary order. 



III. THREE-BODY BETHE-SALPETER EQUATION 
A. Three-body t-matrix and wave functions 

The three-body Bethe-Salpeter equation can be motivated in much the same way as in the 
two-body case. Examining the Feynman diagrams contributing to the three-body t-matrix 
it is possible to separate them into two categories, three-body reducible and irreducible 
diagrams. The irreducible diagrams fall into two classes, those where only two of the three 
particles are interacting and those where all three particles are interacting. The sum of all 
irreducible diagrams can then be written as 

3 3 

V = V° + ^VHG- 1 = V° + ^2V, (14) 

i=i i=i 

where V° is the sum of all three-body irreducible diagrams where all three particles are 
interacting and V % is the sum of all diagrams where only particles j and k (with j ^ k ^ i) 
are interacting. That is V° is a three-body interaction and V 1 is a two-body interaction. 
The three-body t-matrix can then be written as 

T = V - VG° BS T = V - TG° BS V , (15) 

where the three-particle free propagator in terms of the individual propagators reads G° BS = 
—G1G2G3. 

The three-body Bethe-Salpeter six-point Green function is then given as: 

G = G BS ~ G BsVG = G BS ~ SVG° BS 

= g °bs ~ G BS TG BS . (16) 
Obviously, the inverse three-body propagator is 

G- 1 = (CPbsT 1 + V ■ (17) 
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It is convenient to introduce the Faddeev decomposition of the full t-matrix into subam- 
plitudes T ij , 

3 

T=J2 r ' lj > ( 18 ) 

i,j=0 

where the indices i and j characterize the subamplitudes according to the character of the 
last and first interactions; ie., for i ^ (j 7^ 0) , the particles i (j) are not taking part in 
the last (first) interaction, i — (j = 0) means that the last (first) interaction is genuine 
three-particlejirteraction. The subamplitudes T ,J satisfy the following integral equations 
(see also ref. 



r*j = vSij - vg° bs J2 Tkj ■ ( 19 ) 

k=0 

For numerical solution it is more appropriate to re-write these equations in terms of two- 
body and three-body t-matrices which subsume the corresponding interaction kernels j^: 

M l = V i - V i G i BS M i , (20) 
where M l is the two-body scattering matrix and G BS = —iGjGk = iG^ 1 G° BS ; and 

M° = V°- V°G BS M° , (21) 

iterates to all orders the three-body interaction . The subamplitudes for three-body Bethe- 
Salpeter equation for distinguishable particles can be written as 

3 

T iJ = MHGT l 5 i3 -M l G i BS Y J B tk r kj 

k=0 

3 

= M'dij - M'Gls B ik T k i , (22) 

where — 1 — 5^. and Ai l are defined in terms of M l as in eq. ()14|). 

As in the two-body case, the three-body Bethe-Salpeter wave functions can be obtained by 
examining the residues of poles in the interacting three-body propagator. The wave function 
for a three-body bound state with mass M B is related to the residue of the propagator at a 
pole corresponding to P 2 = M B , where P is the total momentum of the three-body system. 
From the resolved form of the propagator, it is clear that this is due to a pole in the t-matrix. 
The three-body bound state vertex function is related to the residue of the t-matrix at this 
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pole so that the t-matrix reads 

\r) (r 
q-i 

Using this in the t-matrix equation, the vertex function can be shown to satisfy the equation 



T = -#r^ + R - < 23) 



|r> = -VG% S |D . (24) 
The interacting three-body propagator can be written as 

G°bs ~ G° BS TIG BS . (25) 



g%s 10 (r| G BS 



Mf - P 2 

The bound-state wave function is defined as 

|*>=G^ 5 |r>, (26) 

which satisfies the wave equation 

Q- 1 \^f>=0. (27) 
The bound state vertex and wave functions are decomposed into their Faddeev components: 

3 

|r > = J2\ ri >> ( 28 ) 

3 

|tf > = ^|^>, (29) 

defined by dynamical equations 

|p > = -v*G^ 5 |r >, (30) 

{Gbs)- 1 \¥ > = -V*|# >, (31) 

where i = 0, 1, 2, 3. 

The scattering states can be obtained from (J16)) . Since has poles at the on-shell 
masses of particles with plane-wave residues, isolating the poles for the right-hand propa- 
gators in (|16J) gives the three-body scattering wave function with outgoing-wave boundary 
conditions 

|tt(+))= (l-G BS T)\p 1 ,s 1 ;p 2 ,s 2 ;p 3 ,s 3 ) , (32) 

while doing so for the left-hand yields the conjugate scattering wave function with incoming- 
wave boundary conditions. 

= (pi,s l ;p 2 ,s 2 ;p 3 ,s 3 \ (1-TG° BS ) , (33) 



where 



\P, s) 



u iPi s ) \p) f° r s pi n 



\p) for spin = 

The Faddeev components of the three-body scattering wave function are given by 



(34) 



{P1P2P3 



Q1Q2Q3 (\-J2 TkiG Bsj 



(35) 



In order to extract the scattering states where two of the particles are bound and the 
other is in the continuum, the propagator eq. (j!6|) is rewritten using (|22|) as 

3 



G = G° BS -G° BS J2^G' 





BS 



j=0 



3 3 



-EX 

t=l j=0 



G BS MHG t G BS 5 l3 - G BS MHG t G BS B lk T kj G BS 



k=0 



(36) 



The last term in this expression contains the two-body scattering matrices which have poles 
corresponding to the two-body bound states and can be represented by 

|p(2)^ ^(2)i| 



Mf - P i 



+ 11 



(37) 



where |r^ 2 ^) is the two-body bound-state vertex function for the pair of interacting particles 
j and k (J 7^ k 7^ i), M 2 is the mass of the two-body bound state and P i is its total 
momentum. This term also contains the free propagator for particle i. By selecting the 
residues of the pole for particle 1 (putting i — 1) and the bound-state of particles 2 and 3, 
the conjugate state two-body scattering state with incoming- wave boundary conditions can 
be identified as 

3 3 



(OK-)) = (T^; Pl , Sl \G\ 



BS 





BS 



= <$ (2)1 ;pi, S i| 

i=0 



k=0 i=0 
3 3 



k=0 i=0 



k^l 



BS 



(38) 
(39) 
(40) 



A similar expression can be obtained for the wave function with outgoing-wave boundary 
conditions. 
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All of the wave functions satisfy the wave equation 

Q- 1 \^f>=0. (41) 

B. Electromagnetic current in Bethe-Salpeter framework 

The three-body Bethe-Salpeter current can be determined by considering all diagrams 
contributing to the seven-point function with six legs corresponding to the three incoming 
and outgoing particles and one photon leg. By separating the diagrams into three-particle 
reducible and irreducible contributions the current operator can be identified as the sum of 
all irreducible seven-point functions. There will be three types of contributions: one-body 
contributions where the photon attaches to only one of the interacting particles, two-body 
contributions where the photon attaches internally to a two-body interaction, and three- 
body contributions where the photon attaches internally to a three-body interaction. 

The one-body currents are of the form 

J* = -JfG^Gj 1 , (42) 

where t — 1, 2, 3, r ^ s ^ t and Jf is the vertex for attaching a photon to particle t, for 
which q^Jj 1 = [e t (q) , G^ 1 ~\ . Each of these currents satisfies the Ward-Takahashi identities 

q l >J*=Mq),(&Bs)- 1 ] ■ ( 43 ) 

As in the case of two-particle system Q|, the two- and three-body currents follow 
from attaching a photon line to all particle lines and into momentum-dependent vertices 
internal to the two- and three-body Bethe-Salpeter kernels. Defining as the two-body 
current as extracted from the two-body problem and j( 3 )^ as the completely connected 
three-body current, the exchange current for the three-body system can be written as 

jt» = [j^Gt\ t = 1,2,3 (r^s^t) 

which satisfies the Ward-Takahashi identities 

q jty _ J M?) + e s (q), V] , t = 1, 2, 3 (r^s^t) 
</ """ x " ~ Mg),V°l, t = 
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FIG. 1: This Feynman diagram represents part a contribution to the seven-point function. The 
particles are label 1 to 3 from top to bottom. The ovals represent two-body scattering matrices. 

where e T (q) = e x {q) + e 2 (q) + e 3 (g). 

Summation of the diagrams to give the seven-point function must be done with consider- 
able care [l(3]. The three-body Bethe-Salpeter equation is based on Feynman perturbation 
theory rather than on time-ordered perturbation theory as in the non-relativistic case. So 
far this difference has been manifest only in the need to introduce inverse propagators into 
some of the operators in order to obtain equations with a similar structure to the non- 
relativistic case and in the fact that the implied integrals are four- dimensional rather than 
three-dimensional. The problem with seven-point function can be seen by considering Fig.QJ 
Here the virtual photon is absorbed on particle 1 and is preceded by two-body interactions 
between particles 1 and 2. Between these interactions is two-body t-matrix for particles 2 
and 3. Since this is a Feynman diagram, there is no ordering implied between the absorbed 
photon and the t-matrix for particles 2 and 3. In summing the interactions this t-matrix 
can either be associated with the sum on the left or with the sum on the right, but not with 
both. As a result, if we associate this interaction with the right hand sum the last interac- 
tion to the right must be between particles two and three while the first interaction on the 
left must be between particles 1 and 2. This requires that there must be a restricted sum 
between three-body scattering matrices which occur on either side of a one-body current. 



The possibility of double- counting this contribution was pointed out in 



111 ]. The seven-point 
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function can then be written as 

0" = " E (g°bs ' G° BS E T ij B jt G% s ) J* ( G BS - G° B5 £ T fc <G' 



o 

BS 



t=l 



i,j=0 



k,l=0 



-s E^K 



(46) 



,i=0 



where the restricted sum is represented by the Bj t associated with the left-hand t-matrix. 
While the exchange current term is nicely factored with complete interacting six-point propa- 
gators on the left and right of the current operator, the one-body current term is not factored 
in this way. This factorization can be accomplished by using the identity 



g- 1 G° BS - G° BS £ T«B jt G° BS ]=!+ V'Gls = 1 + V t iG t 



BS > 



(47) 



i,j=0 



to rewrite the seven-point propagator (|4T)|) as 



g» = -g I J2 ((l + V l G BS ) J*" + JpG" 1 ) + J {3) A g 



-g 



.t=l 

3 



t=l 



t=0 



(48) 



where we have defined the one-body and interaction currents in three-particle space and 
J*^., respectively, by eq. (}4*2j) and by 

Jfi = ^ + iV*JT for t ^ , (49) 

■S = = J {3> , (50) 



since 



(51) 



behaves like a two-body current. Thus, the total Bethe-Salpeter current is identified as 

= E Jt " + E ^ = J (l) " + JL , (52) 



t=l i=0 

and the WT identity for this current can be easily shown to be 



(53) 



This along with the wave equation for the bound and scattering states implies that matrix 
elements of this current between any three-body states is gauge invariant: 



(54) 
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FIG. 2: A diagrammatic representation of the interaction current resulting from the symmetric 
factorization of the seven-point function. The rectangular box represents the two-body Bethe- 
Salpeter kernel and the rectangular box with attached photon represents the corresponding two- 
body electromagnetic current. 

C. Identical particles 

For the case of identical particles it is necessary to explicitly symmetrize the n-point 
functions. In particular, the transition to the corresponding functions for the Gross equation 
is simplified by writing the n-point functions in terms of explicitly symmetrized interactions. 
Any full operator (summed over all three particles) K (e.g. 1Z = V,T, Q, J^, j£ t , J e M ff , . . . ) 
commutes with any particle interchange operator Vif 7Z = VylZVij. Therefore we can write 
for its symmetrized form R 

R = A 3 Tl = TZA 3 = A3KA3 , (55) 
R = A 3 R = RA 3 , (56) 

where the normalized three-body symmetrization operator is 

A 3 = y (1 + CVu + CP13 + (V23 + v A + v 5 ) 
= ^(l + (V l ,)(l + (V tk + (V lk )=A k 2 U k 
= ^ (1 + CP ik + (V jk ) (1 + (V i:i ) = U k A h 2 , (57) 
with £ = 1 (£ = — 1) for bosons (fermions), respectively, and with 

V A = V 13 V 12 = V 23 V 13 = V l2 V 23 , 

V 5 = V 23 V 12 = V l2 V 13 = V l3 V 23 , (58) 
14 



and with A% being the two body symmetrization operator acting in the two-body subspaces 
of particles i and j (where i / j / fc) and the sum of additional exchanges mixing this 
subspace with other channel, which we denote 

A = \ {l + CPij) , (59) 

n fc = i (l + (v ki + (v kj ) . (60) 

Note that A 3 = A 3 A k 2 = A k 2 A 3 for any k. We will exploit this property in the following 
discussion. 

The symmetrized t-matrix reads 

T = A 3 T = V- VG° BS T = V- TG° BS V , (61) 

with V being the symmetrized potential: 

3 

V = A 3 V = A 3 V° + J2 U * V* iG 7 X , ( 62 ) 
i=i 

where V % = A\V % is the symmetrized two particle interaction. The symmetrized six-point 
function can then be written as 

G = A 3 Q = A 3 G° BS - G° BS TG BS = A 3 G° BS - G° BS VG = A 3 G° BS - GVG° BS . (63) 

From eqs. (|61l63j) the symmetrized bound and scattering state wave functions follow as in 
previous section. It is easy to see that for full wave functions it holds 

|* s >= A 3 \^ > . (64) 

One can also define the symmetrized Faddeev components of wave functions: 

> = .4 3 |*° >, (65) 

m> = l(A\^> +av^i* j > +cPikA%\* k > ) , (66) 

where % ^ j ^ k. Making use of the symmetry relations (|B48|) the symmetrized wave 
functions is expressed in terms of only two of these components: 

|* s >= |*2 > +(1 + CP12 + CV 13 )\^l >= |*2 > +311! I** > . (67) 
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The symmetrized seven-point function can be written as 

G» = A a g ,t = -GJ&G, (68) 

and the symmetrized current matrix elements are 

t m =< J» s |* s > = < <°| > +3 < <°| J c M ff |*J > + 

3 < ^1 |*2 > +3 < *; x | (1 + 2CP 12 )|*] > , (69) 

where we have used the symmetry relations for the Faddeev components of wave functions 
and the fact that the current commutes with any replacement operator V%y 

The BS formalism can be conveniently represented in the matrix notation with matrix 
indices corresponding to the Faddeev components. This matrix notation is particularly useful 
for further re-arrangement leading to the spectator formalism. We give a detail description 
of this matrix notation in appendices. 



IV. THE THREE-BODY GROSS EQUATION 
A. T-matrix and wave functions 



The three-body Gross equations can be obtained from the three-body Bethe-Salpeter 
equations by placing two of the intermediate-state particles on their respective positive 
energy mass shells. For amplitudes whose first or last interactions are two-body interactions, 
one of these on-shell particles must be the spectator [5j. For amplitudes whose first of last 
interaction are three-body interactions any pair of on-shell particle can be selected Q]. This 
prescription will work only for three identical particles since it precludes the possibility that 
the on-shell pair can form a bound state. Since the choice of two of the particles to be set 
on shell clearly violates the quantum symmetry of the three-body system, it is necessary to 
symmetrize the interaction kernels to impose the appropriate symmetry of the system, as 
was done in the case of the two-body system [l^ . 

The first step in obtaining the three-body Gross equation is to obtain the three-body 
Bethe-Salpeter equation with symmetrized interactions. This can be done conveniently by 
first expressing the Faddeev components of the various three-body Bethe-Salpeter ampli- 
tudes in terms of four-dimensional matrices and vectors as described in Appendix A. The 
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symmetrization can then also be described in terms of matrices as described in Appendix 
B. Once the symmetrization is accomplished, it can be noted that the various Faddeev am- 
plitudes can all be expressed in terms of only two types of amplitudes allowing the various 
amplitudes to be written in terms of two-by-two matrices where the initial (final) interaction 
is either a three-body interaction, or an interaction between particles 2 and 3. This is the 
2-channel formalism of Appendix B. The three-body Gross equation can now be obtained 
from this two-dimensional matrix form by constraining particles 1 and 2 to the positive- 
energy mass shell. It can be done conveniently using the operator formalism introduced in 

Q. 

As described in ^j, a particle can be placed on shell by picking up only the pole of the 
propagator for this particle when performing energy-loop integrals for the Bethe-Salpeter 
equation. This can be done by carefully examining these integrals for the various contribu- 
tions to a given amplitude or by examining all of the Feynman diagrams associated with each 
amplitude. This can be rather tedious. In Q], we developed an algebraic approach where 
various products of operators are defined such that they reproduce the results obtained by 
choosing poles in the energy integrals. This approach has the additional advantage that, 
by starting from the corresponding Bethe-Salpeter amplitudes it is possible to obtain ex- 
pressions for the Gross equation amplitudes that can be easily approximated to any order 
in various expansion schemes. This was shown in detail for the two-body case in |l4j |. The 
primary step in this approach is to replace the propagators for the on-shell nucleons using 
the identity 

Gi — iQi + AGi , (70) 

where Qi is the operator that places the particle on shell and AG, is the difference between 
the full propagator and its on-shell contribution. Note that AG{ is simply the one-body 
Feynman propagator with the wrong boundary condition. 

Starting with the two-dimensional forms of the Bethe-Salpeter amplitudes, we replace 
the three-body Feynman propagator with the sum of a piece with particles 1 and 2 on shell 
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g\ Qi 0.2 and a remainder Ag\ using 



G» BS = -G 1 G 2 G 3 = -(iQ 1 + AG 1 )(iQ 2 + AG 2 )G 3 

= g 1 2 Q 1 Q 2 + Ag 1 2 , (71) 

g\ = G 3 , (72) 

= (-iQiAG 2 - AG ± iQ 2 - AG X AG 2 ) G 3 . (73) 

Substituting this into the two-dimensional form of the Bethe-Salpeter scattering matrix 
flB84j) gives 

t = v-v8(g 1 2 Q l Q 2 + Ag 1 2 )t, (74) 

where the matrix 8 is defined by (jB76|) . This can in turn be reexpressed as a pair of coupled 
equations as 

t = u - uglQiQ 2 St = u - tglQiQ 2 6u, (75) 
u = v — v Ag^S u = v — u Ag\8 v , (76) 

where u is the quasipotential. 

This system of equations is closed by multiplying from left and right by QiQ 2 which 
leads to 

t\ = ul-u\gl8tl = u\-tlg\8u\, (77) 

where u\ = QiQ 2 uQiQ 2 and t\ = QiQ 2 tQiQ 2 . 

Although (|77|) appears to be much like the three-body Beth-Salpeter equation (jB84|) . it 
is in fact much more complicated due to the quasipotential equation (|75j) . The appearance 
of the matrix 8 in f!76|) means that, in general, the quasipotential matrix is also nondiagonal 



whereas the Bethe-Salpeter kernel v is diagonal. In accordance with (|B61j) we denote 

/ jjOO jjOl \ 

u=\ . (78) 

\U 10 U u J 

The off-diagonal components of u are the result of mixing contributions from the two- and 
three-body Bethe-Salpeter kernels which results in contributions which are also completely 
connected and behave as three-body interactions. The diagonal elements will also be compli- 
cated by 8. The upper left-hand component U 00 will always contain at least one contribution 
from the three-body Bethe-Salpeter kernel in all of its terms, so it will always behave as a 
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(a) 



(b) (c) 
i— +2C — rnhX^v_ 



(d) 



(e) 



— 'V* 1 — — — 

_rW-^ +2C — rWi- + 



2C, 



(0 



(g) 



FIG. 3: Contributions to U up to second order in Bethe-Salpeter kernel. The rectangular boxes 
represent the (symmetrized) two-body Bethe-Salpeter kernel, the crosses represent on-shell propa- 
gators, the open circles represent the off-shell propagators and the filled squares represent factors 
of the inverse one-body propagator. 

three-body interaction even though contributions from the two-body Bethe-Salpeter kernel 
will also occur. 

The lower right-hand contribution U 11 will also be complicated by this mixing. Figure El 
shows diagrams for U 11 associated with the first iteration of (|75|): 

U u ~ l -V\G^ - V^A^l + 2CV l2 )iG^V l . (79) 

The rectangular boxes represent the two-body Bethe-Salpeter kernel for particles two and 
three, the cross represents the on-shell propagator iQi, the open circle represents the off-shell 
propagator AG^ and the square box represents the inverse one-body propagator iG± . These 
diagrams therefore contain some ill-defined products such as QiG^ 1 that can be resolved 
only when the quasipotential is used as part of some larger diagram. The resolution of 
these products is then obtained by using the rules developed in [14( which are repeated in 
Appendix C for convenience. 

Note that in Fig. El particle 1 is disconnected from particles 2 and 3 in diagrams (a-d), 
while the diagrams (e-g) are completely connected. In (|77j). the quasipotential appears only 
with the external legs for particles 1 and 2 constrained to be on shell due to (jClj) and (jC2|) . 
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fa) 



(bl 



(cl 



FIG. 4: Selected diagrams contributing to U 11 at third order in the Bethe-Salpeter kernels. The 
large rectangular box represents the three-body Bethe-Salpeter kernel. 

Using these, the completely constrained quasipotential contains only diagrams (a) and (b) 
which are disconnected and have the form of two-body interactions. As we will see below, 
the quasipotential with only the initial or final legs constrained can contribute to effective 
current operator and diagrams (e) and (f) can appear in that role. Diagrams (c), (d) and 
(g) will never contribute. 

There will be many diagrams that contribute to U 11 at the second iteration of (|76|). 
Examples of such diagrams are shown in Fig. 0] Diagram (a) is disconnected and can 
be categorized as a two-body interaction. Diagrams (b) and (c) are, however, completely 
connected diagrams and are three-body interactions. Therefore, U 11 contains both two- and 
three-body interactions. 

As a practical matter, it may be useful to separate the two- and three-body interactions. 
This can be done by first defining a new, diagonal quasipotential 



w = v — vAg 2 aw = v — wAg, 2 otv 



A 3 W° 







o \w icr 



(80) 



where the diagonal matrix ex. is given by ()B65j) . Thus the upper left-hand component of w is 
an iteration of the three-body Bethe-Salpeter kernel, while the lower right hand component 
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is an iteration of the two-body Bethe-Salpeter kernel: 

W 1 = V 1 - V\-iAG 2 G 3 )W 1 , (81) 

i.e., W 1 is the two-body spectator quasipotential for particles 2 and 3. The component W u = 
^W 1 iGi 1 is identical to the two-body contribution to U 11 . Indeed, u is related to w by 

u = w - wAgl/3u, (82) 

where the matrix (3 is given by (|B78j) . The matrix w can now be written as the sum of two- 
and three-body contributions as 

w = w {2) + w {3) . (83) 
Since is identical to the two-body content of u, u can be written as 

u = w {2) + u {3) . (84) 
Substituting and (jSljl into (JH2J), the equation for the three-body contributions to u is 

u (3) = ^(3) _ ujAgipyjW _ wAglPu® . (85) 

Equation (ffHjl can be resumed to give 

t = m-mQ 1 Q 2 gl/3t, (86) 

where 

m = u — uQiQ 2 g\ot.im . (87) 

Since, in general, u is not diagonal, m is also not diagonal and all of its components contain 
three-body interactions. It is also useful to separate the two- and three-body contributions 
to m. The two-body contribution will just contain iterations of the two-body contribution 
to u. We can then define the two-body scattering matrix as 

= w& - Q x Q 2 glam® . (88) 

Substituting m = + into (|%7j) and using (|HHj) the three-body scattering matrix is 
given by 

m (3) = u (3) _ u (3)g 1 Q 2 ^i am (2) _ u Q 1 Q 2 gla.m { ^ . (89) 
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As in the matrix form of BS equation (cf eqs. (|B86IB87|) ). in the Gross formalism there 
are two propagators, a right-hand one 



9m = - 02*2*02 = 02 1 - 92 U 2 S 92R 



(90) 



with the inverse 

{g\ R )- 1 = {g\)- l l + u\8 , (91) 

and a left-hand one 

al L = oil - 02**202 = gl 1 - 9\ L & u \g\ ( 92 ) 

with the inverse 

(g l 2L y l = (g l 2 r l l + 8ul (93) 

The three-body bound-state vertex function can be obtained using (|77jl and the analytic 
form of the t-matrix at the bound-state pole. The vertex function satisfies the equation 



|T2> = -u\8g\ |^2> = -mkPd IT2) 



(94) 



From the non-linear form of the t-matrix equation in the Gross formalism 

t\ = u\ - t\bg\t\ - t\g\5u\8g\t\ 
we get the normalization condition for vertex function (}94"j) 
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dgl jr U^L 1 



dP 2 " z OP 2 
The bound-state Gross wave function is defined as 



P 2 =M 2 



(95) 



(96) 



W> = & l K> 



(97) 



Using (J54"j) it is straightforward to show that 



0. 



(98) 



The Gross wave functions for the scattering states can be obtained by considering the 
analytic properties of (fl^z,) -1 ) an d then the effective current operator can be determined 
in a similar fashion from the seven-point function (|B105J) and the substitution (j71j) . An 
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alternative approach is to express the Bethe-Salpeter wave functions in terms of the corre- 
sponding Gross wave functions. To do this substitute (fTTj) into the Moller operator giving 



1 - 8tG° BS = 1-St (g 1 2 Q 1 Q 2 + A^ 1 ) = 1 - Hg\Q, x Q 2 - 8tAg\ 
Using (J75|) to replace the t-matrix in the last term gives 

l-StG° BS = 1 - 5tglQ x Q 2 - 6 {u - tglQ x Q 2 8u) Ag\ 
= {l-StglQiQ 2 ) (1-SuAgl) . 

Similarly, the right-hand Moller operator can be rewritten as 



(99) 



(100) 



1 - G° BS tS = (1 - Ag\u8) (1 - g\QiQ 2 tS) 



(101) 



By examining the bound-state pole in the t-matrices, we can find expressions for the Bethe- 
Salpeter bound-state wave functions in term of the corresponding Gross wave functions: 



and 



(il>\ = (il>l\ (1 - duAgl) 



\^) = {1-Aglu8)\^\) 



(102) 



(103) 



The three-body scattering final state can be obtained directly from the two-dimensional 
form of the corresponding Bethe-Salpeter state (|B95|) as 



(-) 



^ (Pi, s%; p 2 , s 2 ; p 3 , s 3 \ A 3 d T [l - 8tG° BS ] 

t (Pi, si; p 2 , s 2 ; p 3 , s 3 \ A 3 d T (l - 8tglQ x Q 2 ) (l - SuAg\) 



where 



(1 - SuAgl) 



2 (Pi^i,p 2 ,s 2 ;p 3 ,s 3 \A 3 d T [1 - 6t\g\] 



is the Gross three-body scattering state which satisfies the wave equation 



i(-) 



{qIl)- 1 



0. 



(104) 



(105) 



(106) 



Finding a similar relation for the two-body final scattering state is a little more compli- 
cated since it contains a two-body bound state which also needs to be re-expressed in terms 
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of the corresponding spectator wave function for two-body bound state as well as rewriting 
a Moller operator of a slightly different form than that examined above. The two-body final 
state in the two-dimensional matrix form is given by (jBlOlj) 



<0 1( -1 = ^ (2) ^i,si|A 3 (o l) (1-/9*3,5) • 



(107) 



The symmetrized two-body spectator bound-state wave function is for particles 2 and 3 is 
defined as 

■- Q 2 A 23 \& 2) ) . (108) 



(2) 
<P2 



and satisfies the equation 



(G 3 



W 



22) 



(2) 
f2 



(109) 



where W 22 — Q2W 1 Q\Q 2 . The corresponding Bethe-Salpeter wave function can be ex- 
pressed in terms of the spectator wave function one by 



A 23 \^) 



-iAG 2 G 3 )W 



(2) 
<P2 



The asymptotic spectator state satisfies the equation 



(2) 



w 



22 



0. 



(110) 



111^ 



and the asymptotic state in ()107j) can be now written in terms of this spectator state as 



,(2). 



($ (2)1 ;Pl,Si|^. 23 = < .p, -:pi..s-i 
From Q x Ag\ = -iQ ± AG 2 G 3 , it follows 



Qi 



l-W (-iQiAG 2 G 3 



Qi 



1 - W 1 (-iQ 1 AG 2 G 3 )\ ( 1 ) = Qi ( 1 ) (1 - awAgl) 



;i!2) 



(113) 



The Bethe-Salpeter scattering state can now be written as 



1 / (2) 

g \<P2 5Pl,Sl 



( 1 ) (1 - awAgl) (1 - f3tG° BS ) (114) 



Using the identity 

(1 - awAgl) (1 - (3tG° BS ) = [l - (l - ocwAgl) f3tQ 1 Q 2 g 1 2 ] (l - 6uAg\) 
this can be rewritten as 



(115) 



!(-)! 



(2) 



( 1 ) [1 - (1 - «^A^) /3tQ 1 Q 2fe 1 ] (1 - SwA^ 1 ) 







i(-) 



(1 - (JuA^ 1 ) 



(116) 
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where 

1 / (2) 

y?2 ;pi,si 



3 

is the Gross two-body scattering state. It satisfies the equation 



( 1 ) [1 - (1 - ctwAgl) ptQxQzgl] (117) 



<t>t ] 



^l)"=0. (118) 



B. Electromagnetic current 



According to previous subsection, the relationship of the Bethe-Salpeter to Gross wave 
functions always involves a factor of (1 — Ag\u8) for initial states and (1 — 8uAg\) for 
final states. Any Bethe-Salpeter current matrix element can now be converted to the cor- 
responding Gross matrix element by defining the effective current for the Gross equation as 

& = Q1Q2 (1 " SuAgl) r cS (1 - Ag\u5) Q 1 Q 2 , (119) 

where j£ s is given by (|B106J) and (|B107jl . 

Contraction of the photon four-momentum with this effective current operator yields 

= Q1Q2 (1 - SuAgl) {e T (q)Sg^ - g L l e T {q)8) (l - Ag\u5) Q l Q 2 

= Q1Q2 [e 3 (q)d (G^l + uS) - (G^l + Su) e 3 (q)8] Q,Q 2 

= e 3 (q)5Q l Q 2 (gl R )- 1 - (gl^Q^Se^q) , (120) 

where the rules of Appendix C have been used to simplify this expression (the proof is similar 

I I 

to that for a two-body system described in detail in Sec. 3.2 of ref. 14]). This, along with 
the wave equations for the Gross wave functions, implies that the current is conserved. 

To bring the expression for the effective current to a form that can be more readily used 
for practical calculation requires use of the non-associative algebra for the on-shell operators 
Qi and the propagators developed in pjj and summarized in Appendix C. The pieces of 
primary concern are those involving the one-body current and can be simplified with the 
following identities: 

QiQ2</ (1)M QiQ2 = QiQ 2 4 i > (121) 

Q^J^Agl = Q ± Q 2 (JfGx + J%G 2 ) , (122) 
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Ag^^Q 1 Q 2 = (dJf + G 2 J%) Q 1 Q 2 (123) 

and 

AglJ^Agl = -(iQ 1 J?G 1 + G 1 J?iQ 1 )AG 2 G 3 -(AG 1 J?AG 1 )G 2 G 3 
-Ad (iQ 2 J%G 2 + G 2 J%iQ 2 ) G 3 - G X {AG 2 J^AG 2 )G 3 
+ (—iQiAG 2 - AG x iQ 2 - AG X AG 2 ) (G 3 J»G 3 ) . (124) 

In these expressions we have used the identity QiJiQi = ()C6j) to give QiJiAGi = QiJiGi 
and AGiJtQi = GiJiQi. The corresponding part of the effective spectator current then 
reads: 

&,ia = Q1Q2 [J^S-du(G l J^ + G 2 J^)6-S(J^G 1 + 4G 2 )ud 

+du(Ag l 2 J^Agl)5u6] QiQ 2 , (125) 



where Ag\J^ tJi Ag\ from the last term is given by eq. 

The interaction current also contains a piece containing the one-body current operator 
that resulted from our symmetric separation of the seven-point function into a current 
operator and six-point propagators: 



1 / 

3 



-1 



Jsym = o I I & J l ( 1 ) ■ (126) 



Contributions from this current can be simplified using the identities: 



QiQ2iV 1 J?Q 1 Q 2 = 0, (127) 



QiQiiV J\Ag\ = -(Q^GJiQ.iV G 2 G 3 



(128) 



Ag 1 2 iV 1 J?Q l Q 2 = -{G l JZQ 1 ){G 2 G z W 1 Q 2 ) 



(129) 



and 



Ag^iV J^Agl = (iQ 1 J^G 1 )(AG 2 G 3 iV'G 2 G 3 ) + (G 1 J^iQ 1 )(G 2 G 3 iV x AG 2 G 3 



+ (AG 1 J?AG 1 )(G 2 G 3 iV 1 G 2 G 3 ) . 



(130) 
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The corresponding effective spectator current is given by: 



3 2 eff.sym ~ Ql Q2 ^ # 



ud 



GiJfi^GaGa 



JfdiV^Ga 



Su 



-\-uS 



A^vVfA^ 1 



5 it 



<5QiQ 2 . 



(131) 



The two-body exchange current part of the interaction current is 



^ = 1(1 )5(oi), 



(132) 



where = iG 1 1 J 23 ^ ■ Contributions from this current can be simplified using the identities 



Q1Q2 J e xQiQ2 — QiQzJ^^Qz ) 



(133) 



Q,Q 2 J^Agl = -iQ l Q 2 J { ^AG 2 G, , 



(134) 



A52 JIxQiQ.2 = -iQ 1 AG 2 G 3 J ( g> 1 Q 2 

and 

A^ 1 J% Ag\ = -Q x AG 2 G 3 J ( 2 f AG 2 G 3 + iAG 1 G 2 G 3 J®' i G 2 G 3 . 
The spectator two-body exchange current is thus 



■72eff,ex2 = Q1Q2 







J 



-T(2)M 



it (5 



















23 



AG 2 G 3 iJ { 23 )lx 



(135) 



(136) 



ijf 3 ^AG 2 G 3 



uS 





Ag^ x Agl 



5 a 



Su 



(137) 



Since the Bethe-Salpeter three-body current is completely connected there is no simpli- 
fication to the contributions from this current: 



J ex 



A,ex 3 = Q1Q2 (1 " SuAgl) 6 \ " ex " \ d (1- Ag\u8) Q 1 Q 2 , (138) 

1 
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The full effective spectator current is the sum of these four contributions 



Clearly, this effective current is quite complicated in its complete form and would corre- 
spond to a l arg e number of Feynman diagrams. For the purpose of comparing these result 
to those of [n]] and it is useful to compare the simplest case where the interaction is 
given by a one-boson-exchange kernel and no three-body contributions are present. 

V. SPECTATOR FORMALISM IN THE ONE-BOSON-EXCHANGE APPROXI- 
MATION 

The expressions presented above for the spectator formalism are exact. However, in any 
practical calculation using either the Bethe-Salpeter or spectator approaches it is necessary 
to approximate the equations by truncating the kernels. This is usually done by expanding 
the kernels in the number of meson exchanges. Indeed, the only solutions that have been 
obtained for the three-body spectator equations at this time are in the one-boson-exchange 
approximation. Care must be taken in truncating the kernel in order to maintain the wave 
equations and Ward identities. In this section, we will simplify the general expressions to 
the case of the one-boson-exchange (OBE) approximation. 

In the OBE approximation the quasipotentials u and w are equal to Bethe-Salpeter 
kernel v and there are no three-body forces or currents. Then the two-channel formalism 
developed above reduces to just one independent channel: there is only one independent 
component of the t-matrix, of the effective current and of the wave functions. 

The general expressions for the wave functions and currents given in the previous section 
can be reduced to this approximation by replacing u and w by 



dnsnnnnznnj: 




,sym 




(139) 




(140) 



This implies that the spectator t-matrix becomes 
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where, using (J77)l with the shorthand notation V — 1 + 2(Vi 2 and v\\ = Q1Q2V 1 Q1Q2, the 



nonzero component of t is given by 



-22 ^ 22 



22 u 3 



22 .-> ' 22 -^~22 Gs'PV- 



11 

22 



(142) 



The propagators (J9*U|) and (jH2J), and their inverses (pTTj) and (J9*3~j) in the subspace of this 
symmetrized Faddeev component read: 



9ir — &3 — 3G 3 T 22 V G 3 , (g 2R 



11 



^ = G 3 - 3G 3 V T\l G3 , (g l 2L Y l = G 3 1 + V V 22 



(143) 
(144) 



The effective current, defined by fjl 251131113711381139^ . can be simplified by setting the 
three-body Bethe-Salpeter current to zero, and keeping only terms involving a single boson 
exchange. This means that in ()125|) . only terms up to first order in u = v can be retained. 
In (|131|) and (|137|) . no terms containing u can be retained since the corresponding Bethe- 
Salpeter currents already contain OBE contributions. Since the three-body current is set 
to zero, J2 eff ex3 = 0- Using the rules from Appendix C, the effective current in the OBE 
approximation is then (still in the two-channel notation): 

/o 



J2eS 



QlQ 2 







\ (V 1 G 2 J% + J^G 2 V 1 - J { 2 f 



Q1Q2. 



(145) 



To reduce it further to the effective operator in the one-channel space of the OBE approxi- 
mation, we have to specify first the corresponding wave functions. 

Since the three-body contributions have been eliminated, the bound-state wave function 
has only one nonzero component, 




where 



(146) 



(147) 



(gin)' 1 W> = o. 

This wave equation is represented diagrammatically by Fig. EK- 

The interaction-dependent part of the wave function (|1(J5|) with the t-matrix (J141)) already 
has only one component, but the free propagation of three particles is partially included in 
the th channel, omitted now in the OBE approximation. Notice that the current ()145j) 
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FIG. 5: Feynman diagrams representing (a) the bound-state wave equation, (b) the three-body 
final-state wave function and (c) the two-body final-state wave function in the one-boson exchange 
approximation . 

always acts on the wave functions both from the left-hand and right-hand side with the 
matrix 6. This actually holds for any properly symmetrized one- and two-body operator, 
the matrix 6 just counts all symmetrization factors. Under the action of 8 the three-body 
scattering wave function in the OBE approximation can be replaced by its one component 
form 



5 4 t 1 - 92*2$) d ^3 \Pl, Si) p 2 , s 2 ; P3, s 3 ) 



1 

V 



A 3 (1 - 9G 3 T 22 ) \p h si; p 2 , s 2 ; p 3 , s 3 ) = 6 



4<- ] 



(1 - 9G 3 T£) A 3 |pi,Si;p2,s 2 ;p3,S3> 



(148) 
(149) 



which follows from identities A3V = 3A 3 and 

1 



4 



-SdA 3 = A* 







1 

V 



(150) 



;isi) 



Also the three-body scattering wave function in the final state multiplied by the matrix 5 
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is just 



i(-) 



5. 



(152) 



Thus in the considered OBE approximation the matrix elements of the current ()145|) should 
be evaluated with the one-component scattering wave functions (|149|) . This wave function 
is diagrammatically represented by Fig. From the relations (|142I144|) the proper wave 
equation 



4 M 



(92 



2LJ 







(153) 



follows, which indicates that the one-component expression ()149|) is the correct three-body 
scattering wave function in the OBE approximation (i.e., not just effectively in the matrix 
elements of operators proportional to symmetrizing matrix 6). This can indeed be verified by 
repeating the whole construction of this paper for the theory without the explicit three-body 
force. 

Some care must be used in truncating the final two-body scattering state (JTT7J)- The 
factor 1 — otwAg^ must also be expanded in order for this wave function to satisfy the wave 
equation in the OBE approximation. In this the factor w must be expanded to one order less 
than the order of expansion. Since we are keeping only OBE terms, w should contain zero 
boson exchanges which means that it should be set to zero. This state then also simplifies 
to the one-component form 



<t>V 



0, 



3 
0. 



>L;pi,si|QiQ 2 [i-6(V 12 t£g 3 ] 



(154) 
(155) 
(156) 



This wave function is represented diagrammatically by Fig. 

Since all wave function have in the OBE approximation only one non-zero component, 
the matrix element of the effective current ()145|) can be with the help of ()151|) rewritten as 
a matrix element of the effective one-component current {J^)\ 2 



m & w) 



m mil 



efl722 



QiQ 2 3 



J£V-V( V 1 G 2 J% 



T7 1 -( 2 ) 
23 



J9G2V — J 



V 



(157) 
Q1Q2, (158) 



where ip an d 4>' are an y of the wave functions discussed above. The Ward identity for this 
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FIG. 6: Feynman diagrams representing the effective current operator for the spectator equation 
in the one-boson-exchange approximation. 



current is 



e*(q),V 
V 



3{[e 3 (g),G 3 1 ] Q 1 Q 2 V + V [e 3 (q),V 



e 2 (q)+e 3 (q),V 



22 



(159) 
(160) 



This, along with the wave equations, shows that the current matrix elements will be gauge 
invariant. 

Figure El provides a diagrammatic representation of all of the contributions to the one- 
boson-exchange Gross current. The matrix elements for various processes can now be pro- 
duced by combining the spectator wave functions with this effective current. Some simpli- 
fication can be obtained by defining off-shell wave functions and scattering matrices. Doing 
this leads to expressions coinciding with the results of Q and Q. 



VI. SUMMARY 



In this paper, we have produced a comprehensive and general description of the three- 
body spectator equations and the corresponding electromagnetic currents starting with the 
three-body Bethe-Salpeter equations and using an operator formalism that has been devel- 
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oped for this purpose. Wave functions for bound state, and two- and three-body states are 
presented along with the effective electromagnetic current operator. The results obtained 
are exact and can be truncated for practical calculation. As an example we have shown 
truncation according to the one-boson-exchange approximation is consistent with previous 
results. 
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APPENDIX A: MATRIX FORM OF THE BS EQUATION: DISTINGUISHABLE 
PARTICLES 

1. Propagators and wave functions 

To make the equation for t-matrix, one can rewrite them by defining the matrix compo- 
nents 

V* = VHG-% = V5 l3 , i,j = 1, 2, 3 (Al) 
v o = v o l = v o 5tQ = v o 5i ^ (A2) 

and the same for matrices M.. In terms of these, the equations for the three-body t-matrix 
reads 

3 3 3 

T ij = V v _ j2 v u g° bs Tkj = - mUg °bs Tkj ' ( A3 ) 

1=0 k=0 1=0 k^l 

where i,j now run from to 3. We will often use the matrix B with components 

B l3 = 1 - 8ij , (A4) 
which allows to re-write the equation for t-matrix in the form: 

3 

T ij = M ij _ M il G° BS B lk T kj . (A5) 

l,k=0 
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The above form suggests a convenient matrix notation in the "channel space" defined by 
respective Faddeev components (T)ij = T % \ i,j = 0, 1, 2, 3. In terms of such matrices 



T = M - MG° BS BT = M- TBG° BS M 

= V - VG° BS (1 + B)T = V- T(l + B)G° BS V , 



(A6) 
(A7) 



with 



M = V - VG° BS M = V - MG° BS V . (A8) 

The complete three-body t-matrix can be related to Faddeev components in the matrix 
form by introducing four- dimensional vector 

/l\ 

1 



D 



1 

V 1 / 



(A9) 



To sum over the channels one has to multiplying the corresponding matrix or vector by the 
vector D T and/or D 

3 



T = 


E 


rpij _ 




i,j=0 




*> = 


D T 




r> = 


D T 





(r| = (t\d. 



(A10) 

(All) 
(A12) 



Substituting this into the six-point function (fTBj) 



Q — G BS - G BS D T TDG° BS , , 



and using 



D D 



[1 + B)D = AD, 



(A13) 

(A14) 
(A15) 



the six-point propagator (fTBj) is re-written as 

1 



A G° BS D T D - G BS D T T- (1 + B) DG% 
-D T [G° BS - G° BS T (1 + B) G° BS ]D 



BS 



(A16) 
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where G° BS = G BS 1. Note that in using ()A15|) in obtaining this result we could have as 
easily used the transpose of this relation which would reverse the order of T and 1 + B. 
We can therefore define two matrix propagators as 

g R = G° BS - G° BS T (1 + B) G° BS (A17) 

and 

g L = G° BS - G° BS (1 + B) TG° BS . (A18) 
An additional useful form of the six-point propagator is obtained by substituting (jA7|) into 

dun 

g R = G° BS - G° BS V (1 + B) G R (A19) 
and similarly from (jA18|) one gets 

Ql = G° BS -Q L (1 + B) VG° BS . (A20) 

The inverses of these propagators are 

Qr = (G^)- 1 + V (1 + B) , (A21) 
Qt = (G BS )- 1 + (l + B)V. (A22) 

By construction D T g RyL D = 4£ and D T g~ R 1 L D = AQ-\ 

The bound and scattering state wave functions can be obtained from consideration of the 
singularities of the t-matrix and propagator. First the three-body t-matrix has simple poles 
at the location of bound states. The matrix form of the t-matrix can then be represented 
as 

ir) (ri 

' ' 1 TZ (A23) 



P 2 - M 2 + ir] 

where \T) is a vector of the bound state vertex functions and 1Z is a matrix of the residual 
parts of the t-matrix. Substituting ()A23j) into (|A7j) and extracting the residues of the 
bound-state poles yields the equation for the vertex function 

|r> = -VG BS (1 + B)\T) , (A24) 
(r| = -(T\(1 + B)G BS V. (A25) 

Defining, the Bethe-Salpeter wave function as |^) = G° BS \T), this can be rewritten as 

[(Gzs)- 1 + V (1 + B)] |*) = g R l |*) = , (A26) 
(*l Wbs)- 1 + (1 + B) V] = <*| g L l = . (A27) 
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The normalization condition for the vertex function can be obtained from the nonlinear form 
of the t-matrix equation 

T = V - TG° BS (1 + B) T - T (1 + B) G° BS VG BS (1 + B) T . (A28) 



Substituting (jA23|) into (jA28|) and isolating the residues of the simple pole yields 

d 



1 = -(r\ — [G° BS (l + B) + G° BS (l + B)V(l + B)G BS ] 
= (T f/ 



r) 



d ^(l + B)-G BS (l + B)^(l + B)G BS 



P 2 =M 2 

|r) 



P 2 =M 2 



Since 



DD T = 1 + B, 



(A29) 



(A30) 



this is equivalent to the normalization condition given in ref. jjj. 
The outgoing scattering wave function (|3l!j) can be re-written as 

(¥<->| = i ( PlP2 p 3 \ D T (1 - (1 + B) TG^) D , 

from which one gets the corresponding vector in the space of Faddeev components: 

1 



(A31) 



( Pl P2 P3 \D T [1-(1 + B)TG° BS ] 



where 



0. 



Similarly, the incoming scattering wave function (}3*2"j) is 



* (+) ) = J [1 " G BS T(1 + B)] B \ PlP2P3 ) 



where 



= 0. 



(A32) 



(A33) 



(A34) 



(A35) 



The two-body scattering state for the three-body system can be obtained from ()38|) as 



($ lH \ = (*® 1 ,p l \(0 1 ) (1 - BTG° BS ) D 



(A36) 



The vector form of this state can then be identified as 



(*K-) | = (*®\p 1 \ ( 1 o)(l-BTG° BS ) . 



(A37) 
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where 



(*K-> | 0-1 = 

The initial two-body scattering wave function is 

/o\ 

1 


W 

where 

g- 1 |#i(+)> = o. 



\* 1 ™)=(l-G» B3 TB) 



|$ (2)1 ,Pi> • 



(A38) 



(A39) 



(A40) 



2. Seven-point vertex function and current operators 



The matrix form of the seven-point propagator can now be defined such that 



g» = D T Q tl D. 



(A41) 



Using this and the properties of D, the matrix form of the seven-point function is given by 
Q» = — L j [1 - G BS (1 + B) T] G BS J^G BS (1 + B) [l - TG BS (1 + B 



+ g °bs (1 +B)[1- TG° BS (1 + B)] J? nt [l - G° BS (1 + B) T] G BS (1 + B) 



^Gl [J {1) » (1 + B) + (1 + B) (1 + B)] Q 



R 



where 



(A42) 



(A43) 



i=l 



and we have introduced the diagonal matrix with components defined by (I49I5U|) : 



int 



diag (J& J& J& J£) , 



%JL = M<?),V] . 



(A44) 
(A45) 



The effective current can then be identified as 



J* = JW" (1 + B) + (1 + B) JL (1 + B) 



(A46) 
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Since J^ t is diagonal, the last term can be identically re-written as (1 + B) J^ t , but the 
form given above is more convenient for discussion below. Contraction of the four-momentum 
transfer with the effective current gives 

g„J& = [erM.GV 1 ] (1 + B) + (1 + B) [e T (q),V] (1 + B) 

= er(q) (1 + B) - Q L X (1 + B) e T {q) . (A47) 

So the current will be conserved. 

APPENDIX B: MATRIX FORM OF THE BS EQUATION: IDENTICAL PAR- 
TICLES 

1. Propagators and wave functions 

To symmetrize the matrix form of the Bethe-Salpeter amplitudes (IA61 IA7|) we consider 
first the diagonal matrices 

Z = diag(Z°, ZHG^ 1 , Z 2 iG^\ Z'^G^ 1 ) , (Bl) 

where Z = V or M. Let us first symmetrize the components of this matrix in respect to the 
exchange of interacting particles defining: 

Z = A ip Z = diagiAsZ^Z^G^^zhG^^iGz 1 ) , (B2) 

where (making use of PjkZ l Pjy. = Z l ) 

T = A\Z l = Z l A\ = A\Z l A\ (B3) 

and the diagonal matrix Ai V is expressed in terms of the symmetrization operators (|57|59|) : 

A ip = diag(A 3 , A\, Aj, A\) . (B4) 

Since 

-A-i.p-Aip <Ai P , (B5) 

we can write 

Z = A, P Z = ZAi P = AipZAip , (B6) 
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and get the symmetrized equations for m-matrices in the form: 



M — V — VG° BS M = V- MG° BS V 



(B7) 



Symmetrizing inside each channel, ie. in respect to interchange of interacting particles, does 
not, of course, account for the total symmetrization, in particular, eg. D T VD ^ V = A3V, 
where V defined in eq. (I62|) and enters as the driving term into equation for symmetrized t- 
matrix (|F)T|) . It is necessary to introduce the matrix A s that symmetrizes between individual 
channels: 

/ 1 \ 

i \CP n |C^i3 

\0 \CVn l(v 23 



/ 



satisfying 



A S D 



( 1 \ 

Hi 

n 2 

where ITj are defined by eq. (|60j). It is easy to see that 

3 

D T A s ZD = A 3 Z = A 3 Z° + u i Z l tG; 1 = Z 



»=i 



as required by eq. (|62j). The fully symmetrized matrix Z is therefore given by 



(B8) 



(B9) 



(BIO) 



Z = A S Z = (A s Aip) Z = AZ 







\z\Gl x \C,V l2 Z 2 iG^ lCVi 3 Z 3 iG^ 



\(V 12 Z iG^ 



IZ^G- 1 



UV 23 Z 3 tG^ 



V 



-{Vy^iG-, 1 \CV 2Z Z\G 2 X \z\Gt ) 



(Bll) 
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where the full symmetrization matrix reads: 

/ A 3 







\ 



3 s\ 2 

\C,V l2 A\ \A\ \C,V 23 A\ 



\ \CV^A\ \QV 23 A, 
( A?, 



1/13 
3^2 



\ 
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\A\CV 12 \A\CV 



13 



\AlCV l2 \A\ \A\C,V 23 
V IAUV 13 IA 3 2 (V 23 \Al ) 
When summed over channels it reduces to the full symmetrization operator A3: 



(B12) 



AD = DAz 



D T A = D T A 3 , 



and with other matrices introduced above it satisfies the identities: 

A 2 = A, 
AB = BA. 

Note that the matrices A s and «-4j P do not commute with the matrix B, although A ■ 
does. From these identities several useful relations follow: 

AZ = ZA=AZA, 

Z = AZ = A S Z = Z A s = A S Z A s = AZ =Z A, 



(B13) 

(B14) 
(B15) 
(B16) 
(B17) 

(B18) 
(B19) 



Now we can symmetrize any general (non-diagonal) matrix in 4-channel space, eg. the three 
particle t-matrix T: 



T = ATA, 

T = A 3 T = D T TD. 

The symmetrized matrix satisfies the identities 

T = AT = TA = ATA 

T y^ipT Ty^ip oAipT ] <As'ip 



(B20) 
(B21) 

(B22) 
(B23) 
(B24) 
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which allows to cast the matrix equation for the symmetrized t-matrix into the form: 

T = M- MG° BS BT = M - TBG° BS M (B25) 



V - VG° BS (1 + B)T = V - T(l + B)G° BS V . 



(B26) 



Since according to (jB19|) Z = A s Z = AZ and matrices A and B commute, one can 
re-write these equations into more convenient equivalent form in terms of diagonal matrices 
M and V: 



TBG° BS M 



MA S - MG BS BT = MA S 
VA S -VG° BS (1 + B)T = VA S - T(l + B)G° BS V , 
VAs - TG° BS (1 + B)T - T(l + B)G° BS VG BS (1 + B)T 



(B27) 
(B28) 
(B29) 



Proceeding exactly the same way as in the previous appendix one gets from these relations 
the dynamical equations and normalization conditions for the symmetrized bound state 
amplitudes: 



|r s ) = a\t) , 

«4.|r s ) A-s A^p |T S ) , 

|r s ) = -VG° BS (1 + B) \Y S ) = -VG% S (1 + B) \T S ) 



dG° BS fl 



dV 

B)-G BS (l + B)-^-(l + B)G BS 



dP 2 y ^ V > Qp2 

The symmetrized bound state wave function is defined by 

\* Bs ) = G BS \T s ) . 



|r s ) . 



P 2 =M 2 



(B30) 
(B31) 
(B32) 
(B33) 

(B34) 



The scattering outgoing states are obtained in a straightforward way acting by the sym- 
metrization matrix A on vectors (jA32IA37|) and using the fact that A commutes with 
(1 + B): 

1 



A 



(p ljSl ;p 2 ,s 2] p 3j s 3 \D T A[l-(l + B)TG BS } , (B35) 



= (^\A=(^ 1 ; Pl ,s 1 \ ( 1 o) A{1-BTG° BS ) 

From eqs. (jA17|A18j) one gets the symmetrized propagators 

G R = G° BS A- G° BS T(1 + B)G BS A, 
G L = G° BS A - G° BS A(1 + B)TG° BS , 



(B36) 

(B37) 
(B38) 
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and from eqs. (|A21IA22|) the inverse propagators 



Gx 1 = AQt = (G£ 5 )-U + V(1 + B), (B39) 
Gl 1 = AQt = (G° 5 )-U+(1 + B)V, (B40) 



which satisfy 



G-^Gr = GrG-J = A, (B41) 
G L l G L = G L G L X = A, (B42) 

and in terms of which the generic equations for the symmetrized wave functions read 

G#|¥.) = 0, (B43) 
(V.lGl 1 = 0. (B44) 

Finally the 7-point function is: 



G^ = AQ 



YqGl [J {1) » (1 + B) + (1 + B) (1 + B)} G R (B45) 



2. Reducing to two independent channels 

The matrix elements of the symmetrized operators and components of the symmetrized 
wave/vertex functions are connected by symmetry relations: 



S 0m = (VijS ™ 


C<mO 


= S m0 (V tJ , | 


*2) 


= CP* 1 


<)l 1 




= (Pij 1 




gir = 


(V jk S ir 


= S ir (V st , | 




= CPjh 


l*i> > 


l r :> 


= CP* 




S ir = 


(V lk S kr 


= S u (V rt , | 




= CVik | 




K) 


= CVik | 





for any m = 0, 1, 2, 3 and i^j^k^0(r^s^t^0). This means that there are only 
four independent components of the symmetrized matrices: S 00 , S m , S 10 , S 11 and only two 
independent components of the symmetrized wave (vertex) functions: ty] (T®,rl). The 
symmetry relations were discussed in Section III of ref. (j3|) (for bound state wave function 
in component form). Let us briefly show how they follow from matrix formalism used in 
this paper. 
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(B49) 



(B50) 



Let us introduce auxiliary matrices: 

Ai ps = diag(£Vij, CP23, CP13, CPn) > 

/ 1 \ 

c u C 12 (Vl2 c 13 CV 13 
C21CP12 c 22 C23CP23 
\0 C31CP13 C32CP23 c 33 / 
Obviously Ai ps Ai P = Ap S for any therefore it holds Ai ps A. = A and thus: 

A PS |* S ) = I*,) , (B51) 
Ai ps S = SAi ps = S , (B52) 

from which the relations (jB46IB47|) follow immediately. The matrix A ss mixes the channels 
1,2,3 arbitrarily. With the help of relations VikPij = VijVjk and A 2 (Vjk = (VjkA 2 = A l 2 it 
is easy to show that 

*A$s<A' *Ass*A.s*A^p *A */4*/4-ss , (B53) 
provided the coefficients c^- satisfy J2i c ij = J2i c ji = 1 f° r an y 3- F° r the matrices A ss with 



so restricted coefficients c^- it then follows 



*A-ss& &*A*ss $ • 



(B54) 
(B55) 



Putting now C\ 2 = c 23 = c 3i = 1 and all others equal to zero, one gets the symmetry relations 
(EH. 



To reduce our 4 channel formalism from the previous section to 2 independent channels, 
taking into account the symmetry relations (jB46HB48|) it is useful to introduce the matrix 

(l \ 
1 

CVi 2 

\ C^13 / 



x 



(B56) 



with the help of which one can write 



I*,) = X\1>) , 

|r s > = x| 7 > , 

S = XsX T , 



(B57) 
(B58) 
(B59) 
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where 




It) 




s 00 s 01 
s w s 11 

The relation ()B61|) holds for all fully symmetrized matrices, e.g., for S = T, M, 
that in particular for operators Z (diagonal for distinguishable particles, e.g., 
follows from eq. f)Bll|) : 

Z = AZ = AJZ = XzX T , 



A*Z° 







\Z iG^ 1 

To invert the relation ()B59|) between S and s we introduce the matrices 

/ 1 \ 
| 
| 
Vo \) 



a 




where X and Y satisfy 



X T YX 



For the symmetrization matrix we get from (jB12|) 
X T Y AYX 



a 



ip 



A 3 
|(1 + CA 3 ) 

To sum over two channels we introduce the 2-dimensional vector T>2'. 



1 ^3 ^ 




{ \A\) 


H 



X T D= \ 


( 1 ^ 




H 




i sn iy 





1 + CVX2 + QV. 



13 



The total symmetrized wave function can now be written as 

|*,) = D T 2 |V) = + (1 + CV12 + (Vis) \*l) , 



44 



and analogously for the vertex function. Similarly, the full t-matrix is obtained by 

T = D T TD = D T XtX T D = D%tD 2 , (B70) 



and analogously for any fully symmetrized matrix. Notice that 

XX T YA = AYXX T = A. (B71) 

Using this relation we easily get 

s = X T Y SYX , S = XsX T , (B72) 
rcxs = X T Y RSYX , RS = XrasX T . (B73) 

Last equation shows that in translating the successive action of two 4-dimensional matrices 
-R and S into 2-channel formalism one has to include the matrix a which simply takes into 
account the existence of the three channels connected by the symmetry relations (JB48J) . 
Notice in particular that the relation A = A A translates into 

a ip = a ip a a ip . (B74) 

The equation for the vertex function and t-matrix contain also the matrices B and 1 + B. 
Although these matrices are not symmetrized (their components do not satisfy (jB46HB48|) . 
they are always multiplied by some fully symmetrized matrix (V, M or V) from which one 
can pull out the symmetrization matrix A. In the transition to the two channel formalism 
we can therefore consider instead the matrices AB and A{1 + B). Since 1 + B = DD T 
and AD = A3D, one gets after simple algebra: 



dip S Clip 1 

(B75) 




8 I 1 , 3 I . (B70) 



3 3(1 + 20P; 



12 



From eqs. (jB67|) and (jB75|) it follows 



X T YBYX = ciip (3 dip, (B77) 

3 
3 6CP12 



,03 

/3 = S - cx = \ \ . (B78) 
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From these relations we obtain for products of the fully symmetrized matrices R or S and 
matrices B and 1 + B the following translation table between 4- and 2-channel formalism 

4channel <-> 4channel 4channel 4channel 

(1 + B)S ^ a ip Ss BS^a vp /3s (B79) 

R(l + B)^r6a w RB^r/3a ip (B80) 

#(1 + B)S ^ r<5s RBS^rj3s (B81) 

3. T-matrix and wave functions in two channels 

Making use of relations from the previous subsection we can now easily re-write the 
equations for t-matrix and wave functions in the formalism with two independent channels. 
Multiplying the equation (jB7|) for the m-matrix with *A S and using (|B62|) gives 

m = v — v G° BS ctm = v — m G° BS a v . (B82) 

For the t-matrix we immediately get (from ()B25HB26|) with the help of relations (|B79IIB8l| : 

t = m-m G° BS f3t = m- tG° BS l3m (B83) 
= v - v G° BS 6t = v - tG° BS Sv (B84) 
= v - tG° BS St - tSG° BS vG° BS St , (B85) 

In a similar way we obtain from eq. (|B37IB38|) the 2-channel propagators: 

g R = X T Y G R YX = g R a ip = (G° BS 1 - G° BS 1 8 G° BS ) a ip , (B86) 
Ql = X T Y G L YX = g L a ip = (G° BS -G° BS St G° BS ) a ip . (B87) 

It is convenient to separate the symmetrization matrix a,i P from the full propagators g R) L 
and define g R> L as above. With the help of (|B84|) one gets for g R ^ 

g R = G° BS 1 — G° BS t S G BS = G BS 1 - G° BS v8g R , (B88) 
g L = G° BS 1 -G° BS dt G BS = G BS 1 -g L 8v G° BS , (B89) 



and for the inverse matrices 



g-^ = l(G BS )- l +v8, (B90) 
g^ = liG^y' + Sv. (B91) 
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For the bound state vertex function we get form (|B32|) : 



| 7 > = -*,<5G° B5 |7) = -m/3G2, 5 | 7 > 



(B92) 



or equivalently 



9R l G BS h)=9 R l W)=0. (B93) 
The corresponding normalization can be obtained either from (jB32|) or from non-linear form 



of the t-matrix equation (JH 

l = (7l 



It) 



P 2 =M 2 

The 2-channel three particle scattering function is defined by 



(-) 



YX 



~(pi,s 1 ;p 2 ,s 2 ;p 3 ,s 3 \ (D T AYX) [l - aa ip dtG BS ] 
j (Pi,si;p2,s 2 ;p3,s 3 \ A 3 d T [l-StG° BS ] 
{Pi,si,p 2 ,s 2 ;p 3 ,s 3 \ A 3 \d T -DltG BS 



where 



and we have used 



d 



T 



1 1 



D T AYX 

A 3 & a ip 
A 3 d T S 

The scattering wave function satisfies 



D T A 3 YX = A 3 d 1 

A 3 1, 

4A 3 DT. 



(-) 



9I 1 



In a similar way: 



4> 



i(-) 



$i(-) 



YX 



($ (2)1 , Pl | (o 1 o) AYX[l-<xa ip (3tG° BS ] 
(^ 1 ;Pl ,s 1 \ (o l)a ip (l-(3tG BS ) , 



(B94) 



(B95) 
(B96) 

(B97) 
(B98) 
(B99) 



(B100) 



(B101) 
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where we have used 



AYX = Xa ip , 

(a 1 o o) x = (o 1 



(B102) 
(B103) 



Also this wave function satisfies: 



> 1( -W = o 



(B104) 



4. The BS electromagnetic current for identical particles 



The symmetrized seven point Green's function (jB45|) is in two-channel notation reduced 



to: 



where 



with 



9" = -^a ip g L [JW"8 + Sj^d] g R a n 



vp 1 



£ s = J^S + 6j? Qt 6 



J mi 



J° M 

o 

u 3 J mt 



From divergences of the one-particle and interaction currents 



o -i 



e T {q),Gl s 



% 8 sL s = [e T (q),SvS] , 



the two-channel form of the total WT identity follows: 



(B105) 
(B106) 

(B107) 



(B108) 
(B109) 



qA = e T(q)Sg R - g L e T (q)8 



(B110) 



APPENDIX C: RULES FOR RESOLVING PRODUCTS OF OPERATORS 



This appendix contains a summary of the rules for resolving products of singular operators 



obtained in 



14] 



Rule 1: 



iQ i G7 1 iQ i -> iQi 



(CI) 
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Rule 2: 

iQi Gr 1 AG, = AG, G'UQ, -> , (C2) 

Rule 3: 

AGiG^AGi - AGi . (C3) 

Note that these identities all refer to products where G" 1 is inserted between factors of 
Z\ = iQi or Zf = AGi, and can be summarized by the compact statement 

Z*G-iZ?^5 u/ Z*. 

For operators O other than Zf, 

Rule 4: 

AG, G" 1 O t = O, G' 1 AG, -> , (C4) 

Rule 5: 

iQ i Gr 1 (9 i = (9 i Gr 1 iQ i -,0. (C5) 

Hence AGj G" 1 — > 1 and iQj G" 1 — > for all operators Oi except Z 1 - . 
The one-body current operator satisfies 

Rule 6: 

Q l J l Q t -> (C6) 

The algebra of these operators is nonassociative so it is necessary to perform these op- 
erations by first applying rules 1-5 to eliminate all occurrences of the inverse one-body 
propagator G" 1 then apply rule 6. 
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